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Abstract 



I A major disagreement between different views about the foundations of quan- 

tum mechanics concerns whether for a theory to be intelligible as a fundamental 
physical theory it must involve a "primitive ontology" (PO), i.e., variables de- 
scribing the distribution of matter in 4-dimensional space-time. In this paper, 
we illustrate the value of having a PO. We do so by focussing on the role that 
Q\ . the PO plays for extracting predictions from a given theory and discuss valid and 

invalid derivations of predictions. To this end, we investigate a number of exam- 
ples based on toy models built from the elements of familiar interpretations of 
quantum theory. 
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1 Introduction 

This paper is based on the view that any viable interpretation of quantum mechanics, or 
more generally any fundamental physical theory, must involve variables describing the 
distribution of matter in space and time. Such variables describe the primitive ontology 
(PO) of the theory [3]. Examples of theories with a PO, examples very relevant to this 
paper, include Bohmian mechanics [21], [121 120] (with the PO given by the particles), 
Schrodinger's first quantum theory [42j H] (with the PO given by the mass or charge 
density), as well as GRWm and GRWf, two versions of the Ghirardi-Rimini- Weber 
(GRW) theory [321 [T5] of spontaneous wave function collapse corresponding to two 
choices of PO: the matter density ontology (GRWm) and the flash ontology (GRWf). We 
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recall the definitions of GRWm and GRWf in Section [2] and that of Bohmian mechanics 
in Section 14.11 

This paper concerns derivations of empirical predictions, i.e., of predictions that can 
be tested empirically. We believe that, for a theory with a PO, a derivation of empirical 
predictions should be based on the PO. To elucidate this statement is the goal of this 
paper. We describe how a valid derivation of empirical predictions should work and 
where it should use the PO and the laws governing it. We do so mainly by means of 
examples. Some of the examples (Section [3]) are valid derivations found in the literature 
concerning GRW or other serious theories, but most of them (Section H]) are novel and 
concern toy theories that we have concocted for the purposes of this paper. Some of 
these toy theories make completely wrong predictions; but that does not preclude them 
from exemplifying valid derivations of predictions (which just happen to disagree with 
empirical findings). For a more systematic analysis of why a fundamental physical 
theory needs to explicitly involve a PO (describing the distribution of matter in space 
and time), see [3j 138] . 

In the following we simply say "predictions" for "empirical predictions," and "em- 
pirical content" for the sum of all empirical predictions of a theory. We do not, in this 
paper, compute any specific predictions for specific experiments. For simplicity, we limit 
our considerations to the non-relativist ic quantum mechanics of N spinless particles; the 
models we describe can easily be modified so as to incorporate spin. 

2 The GRWm and GRWf Theories 

For introductory presentations of the idea behind theories of spontaneous wave function 
collapse, such as GRW theory, see [151 E31 El [301 E]- Detailed introductions to the 
GRWm and GRWf theories have been given recently in [31 HH1 [31] . Here we give only a 
brief description. 

2.1 The GRW Process 

In both GRWm and GRWf the evolution of the wave function follows, instead of the 
Schrodinger equation, a stochastic jump process in Hilbert space, called the GRW pro- 
cess. Consider a quantum system of (what would normally be called) N "particles," 
described by a wave function ip = ip(q 1: ... 7 q N ), q { e M 3 , i = 1,...,N. The GRW 
process behaves as if an "observer" outside the universe made unsharp "quantum mea- 
surements" of the position observable of a randomly selected particle at random times 
Ti,T 2 , . . . that occur with constant rate NX, where A is a new constant of nature of 
order of 10~ 16 s _1 , called the collapse rate per particle. The wave function "collapses" 
at every time T = T^, i.e., it changes discontinuously and randomly as follows. The 
post-collapse wave function ipx+ = hnifX/rV't is obtained from the pre-collapse wave 
function ipT- = lim* *t i>t by multiplication by a Gaussian function, 

ip T+ (q u ...,q N ) = ^g{qi -x) 1/2 ip T -(qi,---,qN), (l) 
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where 

is the 3-dimensional Gaussian function of width a, I is chosen randomly from 1, . . . , N, 
and ^ ^ 

Z = Z(X) =( [ dq x - ■■dq N g(q I - X) |^ T -(<7i, • • -,qN)\ 2 ) (3) 

\Jm.3N J 

is a normalization factor. The width a is another new constant of nature of order of 
10 _7 m, while the center X = X^ is chosen randomly with probability density p(x) = 
Z(x) 2 . We will refer to (Xk,Tk) as the space-time location of the collapse. 

Between the collapses, the wave function evolves according to the Schrodinger equa- 
tion corresponding to the standard Hamiltonian H governing the system, e.g., given, for 
N spinless particles, by 

where m^, k = 1, . . . , N, are the masses of the particles, and V is the potential energy 
function of the system. Due to the stochastic evolution, the wave function ip t at time t 
is random. 

This completes our description of the GRW law for the evolution of the wave function. 
According to GRW theory, the wave function ip of the universe evolves according to this 
stochastic law, starting from the initial time (say, the big bang). As a consequence 
El], a subsystem of the universe (comprising M < N "particles") will have a wave 
function ip of its own that evolves according to the appropriate M-particle version of 
the GRW process during the time interval [£1,^2], provided that ip{t\) = <p(ti) <8> x(ti) 
and that the system is isolated from its environment during that interval. 

Another remark concerns density matrices. It is a standard fact that with every 
probability distribution pi(dip) on the unit sphere 

8(je) = {4> e * M = i} (5) 

of a Hilbert space Jff there is associated a density matrix 



p^ = e^XVI= J (6) 

S(Jf) 

where E means expectation. Since the GRW process is stochastic, xp t is random, and 
with its distribution p t there is associated a density matrix p t = p^ t . It turns out 
that pt evolves according to an autonomous equation (i.e., one that depends on p t only 
through p t ), the "master equation." This equation is of a type known as a Lindblad 
equation, or quantum dynamical semigroup, and tends to evolve pure states into mixed 
states. Although the mathematical details play no role in this paper, we give the master 
equation for the sake of completeness: 

*£t = -j i [H,p t ] + \'E I d 3 xgt /2 x ptgl /2 x -N\ Pt . (7) 

1 j 1 *J 
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Here, g\j 2 is the multiplication operator by the function g(q k — x) 1//2 with g the 3- 
dimensional Gaussian function as in ([2]). 

We now turn to the primitive ontology (PO). In the subsections below we present two 
versions of the GRW theory, based on two different choices of the PO, namely the matter 
density ontology (GRWm in Section l2~2l) and the flash ontology (GRWf in Section l2~3l) . 

2.2 GRWm 

GRWm postulates that, at every time t, matter is continuously distributed in space with 
density function m(x,t) for every location x G M 3 , given by 

N f 2 

m(x,t) =^mi / dqi---dq N 5 3 (qi-x)\ip t (q h ...,q N )\ (8) 

N f 2 

= ^2 m i / dq 1 ---dq i _ 1 dq i+1 ---dq N \il) t (q l} ...,q i _ l} x } q i+1} ... } q N )\ 2 . (9) 

1=1 K 3(JV-1) 

In words, one starts with the |^| 2 -distribution in configuration space R 3iV , then obtains 
the marginal distribution of the z-th degree of freedom qi e R 3 by integrating out all 
other variables qj, j ^ i, multiplies by the mass associated with q iy and sums over i. 
Alternatively, ([8]) can be rewritten as 

m(x,t) = (if> t \M(x)\i/) t ) (10) 

with M(x) = Y2i m i 8 3 (Qi—x) the mass density operator, defined in terms of the position 
operators Q^(<?i, • • • , = ft ^(li, ■ ■ ■ , Qn)- 



2.3 GRWf 

According to GRWf, the PO is given by "events" in space-time called flashes, mathe- 
matically described by points in space-time. What this means is that in GRWf matter 
is neither made of particles following world lines, nor of a continuous distribution of 
matter such as in GRWm, but rather of discrete points in space-time, in fact finitely 
many points in every bounded space-time region. 

In the GRWf theory, the space-time locations of the flashes can be read off from the 
history of the wave function: every flash corresponds to one of the spontaneous collapses 
of the wave function, and its space-time location is just the space-time location of that 
collapse. The flashes form the set 

F = {(A 1 ,T 1 ),...,(A fc ,T fc ),...} (11) 

(with Ti < T2 < . . .). Alternatively, we may postulate that flashes can be of N different 
types ("colors"), corresponding to the mathematical description 

F = {(X 1 ,T 1 ,I 1 ),...,(X k ,T k ,I k ),...}, (12) 
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with 1^ the number of the particle affected by the fc-th collapse. 

Note that if the number iV of degrees of freedom in the wave function is large, as in 
the case of a macroscopic object, the number of flashes is also large (if A = 10~ 15 s _1 and 
N = 10 23 , we obtain 10 8 flashes per second). Therefore, for a reasonable choice of the 
parameters of the GRWf theory, a cubic centimeter of solid matter contains more than 
10 8 flashes per second. That is to say that large numbers of flashes can form macroscopic 
shapes, such as tables and chairs. That is how we find an image of our world in GRWf. 

We should remark that the word "particle" can be misleading. According to GRWf, 
there are no particles in the world, just flashes and a wave function. According to 
GRWm, there are no particles, just continuously distributed matter and a wave function. 
The word "particle" should thus not be taken literally (just like, e.g., the word "sunrise"); 
we use it only because it is common terminology in quantum mechanics. 

3 Predictions and Primitive Ontology 

In the PO view, a satisfactory theory should have a PO. In this view, the PO also 
permits the derivation of predictions. For example, if we want to derive that in a 
certain experiment the pointer of the apparatus will end up pointing to the value z 
with a certain probability then, according to the PO view, we need to derive that the 
configuration of the PO will be such that the matter of the pointer is in a configuration 
corresponding to the pointer pointing to z. In contrast, it would not be appropriate in 
this view to merely show that the wave function lies (approximately) in a subspace of 
Hilbert space corresponding to the pointer pointing to z. 

For example, a number of empirical predictions of GRW theory have been derived 
in J32J SSI ESI ESI [T]. It was found that the predictions deviate from those of quantum 
mechanics but only so slightly that no experimental test has been possible so far [TJ [29] . 
However, the logical clarity of the derivations in (321 EH ESI ESI EE] leaves something to 
be desired, as they do not refer to the PO but limit themselves to analyses of the wave 
function. While we do not dispute that the claimed predictions are indeed predictions 
of GRWm and GRWf, we do see a gap in the derivation. The situation is similar to 
that of a calculation that yields the correct result but is not mathematically rigorous. 
Here, the problem is not one of mathematical rigor but of clarity — philosophical, onto- 
logical, conceptual, and physical clarity. We will describe in this chapter, particularly 
in Section I3.4[ how to close this gap. 

3.1 Calibration Functions 

Given an experiment S ', its outcome Z is a function of the (configuration of the) PO, 

Z = ((PO). (13) 

That is, in GRWf Z = (cRWf(F), an d in GRWm Z = CGRWm(^)- Similarly, the PO 
plays a key role for the claim of empirical equivalence between two theories (for which 
there will be several examples in this paper). 
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Let us elaborate a bit on what the function £ should look like. In GRWm, it is 
natural that Z should be a functional of m(-,t), the distribution of matter at the time 
t when the experiment is completed. What the ( function does is essentially to read off 
the outcome from the display of the apparatus. For example, if the outcome is displayed 
by means of the position of a needle on a scale, ( should read off from m the position of 
the needle, and may concretely be the following: Suppose the experiment is so arranged 
that the region R C M 3 , for simplicity a cuboid R = [a 1; b\] x [a 2 , b 2 ] X [a 3 , b 3 ], contains 
no other matter than the needle, and suppose the scale is along the xi-axis between a\ 
and 61. Then the mean x\ coordinate of the matter distribution at time t inside R is 
given by 

Ldx 1 dx 2 dx 3 x 1 rn(xi,X2,x 3 ,t) 

v^ 1 = r a a a 1 7V~ ' ( 14 ) 

J R ax 1 ax 2 ax 3 m[xi, x 2 , x 3 , t) 

and a typical choice of calibration function would be 

C(m) = z + a{xi) (15) 

with suitable proportionality constant a. Note that ( is a functional of m(-,t). If we 
require that Z be a discrete variable, we may replace (fl~5l) with a suitable step function, 
such as 

C(m) = z + [a{xi)\ , (16) 

where [z] denotes the nearest integer to the real number z. 

In GRWf, Z must depend on the history in an entire time interval, say [t, t + At] 
with t the time when the experiment is completed and At, say, a millisecond. In the 
example of the needle pointing to the outcome, ( needs to read off the position of the 
needle from the flashes and rescale it appropriately. As a concrete example, suppose 
again the region R contains no other matter than the needle, and suppose the scale is 
along the xi-axis between a\ and b\. Then the mean X\ coordinate of the flashes during 
[t, t + At] in R is given by 

Y.k U<T k <t+At lx k eR 

where lc is 1 when the condition C is satisfied and otherwise. A typical choice of 
calibration function would be given in terms of this (x\) by (|T5|) or (|16p . 



3.2 Taking the PO Seriously 

As a consequence of the view that the PO represents matter, we are forced to take the 
PO seriously. 

For example, the transformation behavior of the PO under a Lorentz (or Galilean) 
transformation is constrained by its geometrical nature, and that of the non-primitive 
ontology (i.e., the remaining part of the ontology) is constrained by its relation to 
the PO; see section 4.2 of [5] for elaboration. If the PO consists of flashes, then the 
flashes have to transform like space-time points under Lorentz transformations. If the 
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PO consists of particle world lines, then they have to transform like world lines, i.e., by 
transforming every space-time point on the world line. Put differently, a world line must 
be an unambiguous set of space-time points, independently of the choice of coordinates. 

As another example for what it means to take the PO seriously, we look at a difficulty 
that arose in [6] with an over-simplified discussion of how to derive in GRWm that, in 
a suitable situation, a pointer is pointing to z. The coordinates x\, . . . , x^ € K 3 of the 
configuration space M. 3N of a pointer consisting of N "particles" were decomposed into 
the center of mass x cia = m^Xi/ rrij G M 3 and relative coordinates r ll . . . , r 3N _ 3 , 
and it was derived that in this decomposition the wave function, in a suitable situation, is 
approximately a product ip C m(x cm ) ip te i(ri, . . . , r 3 Ar_ 3 ), where the first factor, the center- 
of-mass wave function ip cm , is a very narrow wave packet. Therefore, it was suggested, 
the matter density associated with the center-of-mass, m cia (x,t) = \tp cm (x,t)\ 2 , is very 
narrow, too (thus making the position of the pointer sufficiently precisely defined). The 
difficulty with this argument is that m cm is not the right quantity to look at: It is not 
the PO, not the matter, not real; it is just a mathematical quantity. What counts is 
whether m as defined in (jSJ) is concentrated in the right location, and there is no simple 
relation between m and m cm . For example, the width of m cm is, in realistic examples, 
10 -13 m whereas that of m is 10~ 3 m (= the width of the pointer). So this argument 
did not take the PO seriously. Rather, it still treated if) in a more or less conventional 
way as providing probabilities for configurations of particles, as that is the situation in 
which the extreme narrowness of ip cm can be relevant. 

The issue of taking the PO seriously also arises in the context of studying the limita- 
tions on knowledge in GRWm and GRWf, a topic we plan to discuss in detail in a future 
work [21] (and that we have outlined in [M]): Inhabitants of a GRWf or GRWm world 
cannot measure the times and locations of the collapses, although these values are well 
defined according to these theories. That is, some things that are real cannot always be 
measured with arbitrary accuracy. But if something cannot be measured, one may be 
tempted to not take it seriously. So, if there are limitations to measuring the variables 
representing the PO, one may be tempted not to take the PO seriously. Needless to 
say, we recommend resisting this temptation. The conclusion "unobservable, therefore 
unreal" is, of course, not a good one. 

Another remark concerns the fine difference between the names "matter density" 
and "mass density" for the m function. The fact that its definition (jSJ) involves the 
quantity rrii usually called "the mass of particle i" suggests the name "mass density" for 
m, which, as we shall argue, can be misleading. We prefer the name "matter density" 
because it reflects what we take to be the fundamental meaning of the m function. Let 
us explain. 

Suppose we had postulated, instead of (jSJ, the following formula for the m function: 

N f 2 

m(x, t) = ^ ej I dq 1 ---dq N 5(q i -x)\i) t (qi,...,qN)\ ■ (18) 

i— 1 

This is the same equation as (|S]), but with m i; the "mass of particle i" replaced by 
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e», the "charge of particle £." In this case, it would evidently not be appropriate any 
more to call the m function the mass density; rather "charge density" would seem more 
appropriate. At the same time, we may postulate that the m function still represents 
the matter density — i.e., we may postulate that matter is distributed continuously with 
density m — and regard (Tl8l) as merely a different law for how matter is distributed!]] In 
that case, the meaning of the m function would not have changed in any way — we would 
only have changed the law governing it. That is why this meaning is better conveyed 
by the name "matter density." Moreover, the matter that we postulate in GRWm and 
whose density is given by the m function does not ipso facto have any such properties 
as mass or charge; it can only assume various levels of density. For example, the m 
function is not a source of an electromagnetic field. 

In addition, the terms "mass density" and "charge density" may easily suggest that 
both quantities physically coexist: Is it not natural to expect that an extended object 
possesses both a mass density and a charge density (different from each other)? In fact, 
it seems completely coherent (though perhaps not desirable) to postulate that the PO 
involves both a mass density and a charge density (different from each other) . If, in con- 
trast, we postulate that the PO consists of only one density — the matter density — and 
choose either (jSj) or (1181) to determine it, then only one of the two functions represents 
something real (more precisely, represents the PO), whereas the other is a pure mathe- 
matical fiction. In other words, in GRWm with (JSJ) as the law of m, the formula on the 
right hand side of (fTBj) lacks any physical significance, just as it does in GRWf. 

3.3 Examples From the Literature 

We now look at earlier arguments based on the connection we are discussing between 
empirical predictions and PO. These are also examples of valid derivations of empirical 
predictions that do not not suffer from the gap we complained about (in the second 
paragraph of this chapter). 

1. The first example is a proof of no-signaling (i.e., the impossibility of transmitting 
messages faster than light between two distant observers, each acting on and ob- 
serving one of two entangled quantum objects) in GRWf due to Bell [15]. The 
proof shows that for two non-interacting systems (here, the system can be taken 
to be one object together with a nearby apparatus), the marginal distribution of 
the flashes pertaining to system 1 depends on the entangled wave function only 
through its reduced density matrix (with system 2 traced out); nor does it depend 
on the Hamiltonian of system 2. Thus, it does not depend upon any message 
that observer 2 may wish to transmit. Now, since outcomes of experiments are 

1 It may be a difficulty with this postulate that (fT8|) can be negative, while it is normal to think of 
matter density as non-negative. This difficulty might be addressed by adding a sufficiently large positive 
constant on the right hand side, or perhaps just by insisting that matter density can be negative, for 
example representing some kind of "anti-matter" (not necessarily related to the usual meaning of this 
word). For the purpose of this section, readers may ignore this difficulty. 
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functions of the flashes, the outcome that observer 1 sees cannot depend on the 
message observer 2 may have wished to transmit, qed. In Bell's words [To] . 

Events in one system, considered separately, allow no inference about 
[. . . ] external fields at work in the other, . . . nor even about the very 
existence of the other system. There are no "messages" in one system 
from the other. The inexplicable correlations of quantum mechanics do 
not give rise to signalling between noninteracting systems. 

A similar proof was given in [4"3T H5] for a relativistic version of GRWf. 

2. The second example is provided by derivations of empirical predictions from Bohmian 
mechanics, which can be found in many papers, e.g., [2JJ Q21 ESI 12] • The PO of 
Bohmian mechanics consists of particles and their trajectories, and the outcomes 
of experiments are read off from the particle trajectories, and not directly from 
the wave function. The empirical content of Bohmian mechanics agrees exactly 
with that of the quantum formalism "whenever the latter is unambiguous." 

3. The third example is taken from the study of superselection rules by Colin et 
al. [22]. For superselection rules it is crucial that some self-adjoint operators are 
not observables. In the terminology of Colin et al., a "weak superselection rule" 
means that no experiment can distinguish between a superposition of vectors from 
different superselection sectors in Hilbert space and a suitable mixture thereof. 
This occurs only if all observables commute with all projections to superselection 
sectors, because otherwise a quantum measurement of such an observable would 
distinguish between them. Colin et al. proved certain weak superselection rules for 
Bohmian mechanics (and Bohmian versions of quantum field theory) as well as for 
GRWm and GRWf, and the proof uses that the outcome of any experiment is a 
function of the PO (and not of the wave function). Indeed, while the PO does not 
permit us to distinguish between a superposition and a mixture of contributions 
from different sectors, the wave function would trivially permit this. 

In more detail, Colin et al. first showed that (in certain situations) the distribution 
of the flashes in GRWf does not distinguish between a superposition and a suitable 
mixture ("strong superselection"), and concluded from this that the distribution 
of outcomes of experiments does not distinguish between them either. Due to the 
empirical equivalence between GRWf and GRWm, the same experiment does not 
distinguish between them in GRWm. To sum up, the proof that some operators 
are not observables used that results of observations must be read off from the 
PO, and then proceeded with a suitable analysis of the PO. 

4. Colin and Struyve discussed in [23J whether their Dirac sea model is empirically 
equivalent with orthodox quantum field theory. They suggested an unusual PO 
containing a huge (or even infinite) number of particles with trajectories, one for 
every electron in the Dirac sea. For the theory to make the right predictions, we 
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need that macroscopic facts such as outcomes of experiments can be read off from 
this PO — and that is what Colin and Struyve analyzed. 



5. In [M] we have shown (and outlined before in [3]) that GRWm and GRWf are 
empirically equivalent. This is a non-trivial statement if and only if the PO is 
taken seriously. 

6. Feldmann and Tumulka (29] considered other values for a and A than suggested 
originally [32] (cr = 1CT 7 m, A = 10~ 16 s _1 ) and evaluated for which points in 
the cxA-plane the GRWm and GRWf theories are empirically adequate, and for 
which they are philosophically satisfactory. Their criterion for being philosophi- 
cally satisfactory is that the PO looks macroscopically like what humans normally 
imagine reality to be like. As a consequence, for determining the philosophically 
satisfactory region they need to pay attention to the behavior of the PO. 

7. In [31] we have derived a measurement formalism for GRWm and GRWf, i.e., a set 
of rules for computing the predictions of GRWm and GRWf for any experiment. 
We call it the GRW formalism; it is spelled out in Section 13.51 below. The GRW 
formalism is analogous to the quantum formalism, i.e., the usual set of rules for 
making predictions in quantum mechanics. Both the quantum and the GRW 
formalism can be formulated as three rules, the first specifying how an isolated 
system evolves, the second specifying the probability distribution of the outcome 
of a given experiment, and the third specifying the change on the system's state 
depending of the outcome obtained. That is, the third rule is a "collapse rule," 
concerning, however, not spontaneous collapses but the collapse induced by the 
intervention of a macroscopic apparatus. Put differently, the GRW formalism is 
what the quantum formalism needs to be replaced with if we pay attention to the 
small deviations from the quantum formalism predicted by GRWm and GRWf. 

The derivation of the GRW formalism in [31] takes the PO seriously. We discuss 
it to some extent in the next two subsections. 



3.4 The Main Theorem About Operators in The GRW For- 
malism 

Before talking about the full GRW formalism, we begin with a statement that is part 
of the GRW formalism, the main theorem about operators: With every experiment S on 
a system and every possible outcome z of <§ there is associated an operator P z acting 
on the Hilbert space J^ SYS of the system. When <§ is performed on a system with wave 
function if), the outcome Z is random with probability distribution 

F(Z = z) = (iP\P z \iP). (19) 

This statement is actually true in both the quantum formalism and the GRW formalism; 
however, to the same experiment, the GRW formalism may assign operators P 2 GRW 
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different from the P^ n assigned by the quantum formalism. Formulas for P® u and P^ GRW 
are given in [31]. (We assume that every experiment has a discrete set of outcomes z.) 

The operator P z in f TIT?]) may be a projection, but in general it is merely a positive 
operator (even in the quantum formalism). If the P z are projections for all z, and if 
the z are real numbers, then the family {P z } corresponds to the self-adjoint operator 
A = J2 z zP z , usually (but misleadingly) called "the observable measured" by S. The z 
are the eigenvalues of A, and the P z the projections to the eigenspaces. When the P z 
are not projections, the family {P z } forms a positive-operator-valued measure (POVM) 
(meaning, in this discrete case, that the positive operators P z are such that J2 Z = 
the identity operator), a concept well known in quantum information theory. 

As an example of a valid derivation of predictions that takes the PO seriously, we 
now outline the derivation of the main theorem about operators from GRWf, following 
[31]. The key point is that we take the outcome Z of the experiment to be a function 
of the pattern F = {(Xi, T 1; 7i), . . . , (X k , T k , I k ), . . .} of flashes, 

Z = ((F) . (20) 

(It would even be realistic to assume that Z depends only on the flashes of the apparatus 
during a short time interval appropriate for reading off the outcome, and not those of 
the system or of the rest of the world, and not those at earlier or later times; but this 
restriction is not needed here.) 

It is a known fact [H] that the joint distribution of all flashes after time t, conditional 
on all flashes up to t, depends quadratically on \l/ t , the wave function of the universe at 
time t. Explicitly, it is given by 

¥(F >t e S\F< t ) = (%\G t (S)\*t) (21) 

with S any set of flash histories after t and GA-) a suitable POVM on the space of all 
flash historiesH 

Let t be the time at which the experiment begins. Consider splitting the universe 
into a system (the object of the experiment) and its environment (the rest of the world, 
including all relevant apparatuses of the experiment), corresponding to a splitting of the 
Hilbert space into = J4? sys <S> if cnv . We assume independence between the system 
and the environment immediately before t, so that 

^t = V®0- (22) 

Here we regard <fi as fixed, while if), the initial state of the system upon which the 
experiment is performed, is allowed to vary in the system Hilbert space J^ ys . We may 
think of as part of the characterization of the experiment, although in practice a 

2 On a continuous space ft, a POVM associates by definition a positive operator G(S) with every 
(measurable) subset S C O in such a way that G(Q) = I and G(Si U £2) = G(Si) + G^) when 
Si H 52 = (and likewise for countable families of pairwise disjoint sets). For discrete f2, G(S) can be 
expressed as ^Z Z&S P Z with P z = G({z}). 
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repetition of the experiment will not begin with exactly the same wave function of the 
apparatus. 

Therefore, the distribution of the random outcome Z is given by 

F(Z = z)= F(F E C\z)) = (y t \G t (C\z)) \9 t ) = (^ GRW |^> , (23) 

where the first scalar product is taken in the Hilbert space of the universe and the second 
in ^ ys , and Pf RW is the POVM given by 

P, GRW = (0|G / (C- 1 (z))|0), (24) 

where the scalar product is a partial scalar product in the Hilbert space of the environ- 
ment. Thus, for every experiment in GRWf, the distribution of outcomes is given by a 
POVM P^ , which is what we wanted to show. 



3.5 The GRW Formalism 

The full GRW formalism is best formulated in terms of density matricesJl Perhaps 
the most remarkable fact about the GRW formalism is that its abstract structure is 
identical to that of the quantum formalism. It consists of three rules, the first saying 
how the density matrix p t of a system evolves when the system is isolated (or, when "the 
observer" is not "taking measurements"); the second (a form of the main theorem about 
operators) saying that with every experiment (or "measurement") there are associated 
operators, and that the probability for obtaining a particular outcome is given by a 
trace formula involving p t ; and the third saying how p t should be changed after the 
experiment, depending on the outcome. Here are the three rules, in both the quantum 
and the GRW version: 

(i) A system isolated from its environment has at every time t a density matrix p t 
which evolves in the quantum formalism according to the unitary (Schrodinger) 
evolution, which for a density matrix reads 

^ = -i[H, Pt } (25) 

with H the system's Hamiltonian, and in the GRW formalism according to the 
master equation (j7]). 

(ii) With every experiment § on a system and every possible outcome z of § there is 
associated a positive operator P z acting on the Hilbert space Jif S ys of the system. 
When $ is performed on a system that has density matrix p at the beginning of 
£*, the outcome Z is random with probability distribution 

F(Z = z) =tr(pP z ) . (26) 

3 That is because we allow the system under study, system 1, to be entangled with another system, 
system 2, which does not interact with either system 1 or the apparatus of the experiment; thus, 
p = tY2\ip) {tp\, where tr2 denotes the partial trace and ip the joint wave function of systems 1 and 
2. Since i\j is usually not precisely known, it is also often convenient to take it to be random; then 
p = Etr 2 |^><# 
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(iii) When the outcome Z — z, the density matrix p of the system gets replaced by 

d = -^BL (27) 

with ^ a (completely positive) linear operation on density matrices. 

In the quantum formalism for an ideal quantum measurement, the P z are projec- 
tions, and 

%(p) = P zP P z . (28) 

Again, the operations ^ 2 ^ u provided by quantum formalism and ^ 2 GRW by the 
GRW formalism may differ. Formulas for c €'^ a and <^p RW are given in 



4 A Set of Examples 

A new and useful perspective on GRW theories arises from contrasting them with other 
theories, even unreasonable ones. For this purpose we develop in this section a set of 
example theories which we obtain by combining elements of the known theories GRWm, 
GRWf, and Bohmian mechanics in new, sometimes playful, ways. Some of the theories 
obtained in this way make completely wrong predictions but are instructive nonetheless 
since they illustrate the way in which predictions follow from a theory. Others make 
predictions in agreement with known empirical facts, yet nobody would seriously propose 
them as fundamental physical theories; still, they allow for illuminating comparisons 
with GRWm and GRWf. 

Here is a "theory construction kit." Choose one of the three primitive ontologies of 
GRWm, GRWf, and Bohmian mechanics: continuous matter density, flashes, or particles 
with trajectories. Choose whether the laws governing the PO should involve a wave 
function tp or a density matrix p. Then choose an evolution law for ip or p, e.g., the 
unitary Schrodinger evolution, the stochastic GRW evolution for if>, or the GRW master 
equation ([7]) for p. Then consider simple laws for how ip (or p, respectively) may govern 
the PO. In this way we arrive at about ten new theories. 



4.1 Bohmian mechanics 

To review the elements used in the theory construction kit, let us briefly recall the laws 
of Bohmian mechanics. Bohmian mechanics is a (non-relativistic) theory of particles in 
motion. The motion of a system of N particles is provided by their world lines 1 1— > Qi(t), 
i = 1, . . . , N, where Qi(t) denotes the position in M 3 of the i-th particle at time t. These 
world lines are determined by Bohm's law of motion J2TJ [121 120] , 

-7T = v i(Qi> • • > Qn) = — Im [Qi . . • , Q N ), (29) 

at rrii ip*ip 

where the wave function if> evolves according to Schrodinger's equation 

ifbjL = Hip , (30) 
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with H the usual nonrelativistic Schrddinger Hamiltonian; for spinless particles it is of 
the form 

An important probability distribution in Bohmian mechanics is the quantum equi- 
librium distribution 

p^q) = mq)\ 2 - (31) 

(While the distribution density is usually denoted p, we write p here in order to reserve 
the letter p for density matrices.) As a consequence of Bohm's law of motion (129 p and 
Schrodinger's equation ( 15U|) . \ip\ 2 is equivariant. This means that if the configuration 
Q{t) = (Qi(t), . . . , Qjv(i)) of a system is random with distribution \ifj t \ 2 at some time 
t, then this will be true also for any other time t. Because of equivariance, Bohmian 
mechanics reproduces the predictions of the quantum formalism for typical initial con- 
ditions of the universe, as discussed in detail in [27J EE] • 

4.2 Bohmian Trajectories and GRW Collapses 

We begin with the particle ontology, with the particle trajectories governed by Bohm's 
law of motion. We consider several ways of combining this with the GRW evolution of 
the wave function or a similar one. 

4.2.1 Bohm's Law and GRW's Law 

First, suppose simply that the particles move according to the usual Bohmian law of mo- 
tion ( )29|) . but that i/j = if) t is the GRW wave function, so that the GRW process replaces 
the unitary evolution. In this theory, which we denote GRWpi, the \ip\ 2 distribution for 
the configuration is not equivariant. 

A world governed by this theory GRWpi has little resemblance with our world: 
It behaves in a very unstable way. For example, a system with the wave function 
of Schrodinger's cat has, even before the collapse into either |dead) or | alive), a con- 
figuration of either a dead or a live cat — but the collapse need not agree with that 
configuration. The cat could be alive before the collapse (i.e., its particle configuration 
is that of a live cat), and still the collapse could reduce ip to a state vector close to 
|dead) — with however no immediate change in the particle configuration of the cat. 

At this point, the reader may feel unsure whether to conclude that just after collapse 
the cat is really alive or that the cat is really dead. That is, partly, what makes this 
theory philosophically useful, despite its empirical inadequacy: It nicely illustrates the 
role of the PO. Taking the PO seriously, we must conclude that the cat is really alive; 
after all, the PO in this theory consists of particles, and the particle configuration is one 
of a live cat. This illustrates that the mere fact that the wave function is one of a dead 
cat does not, in and of itself, mean that there is a dead cat. 

From that fateful collapse onwards, the configuration is guided, in the sense of 
Bohm's law of motion ( 12"9"|) . by that tiny part of i/j that remains of | alive) after the 
collapse, and who knows what happens then. For sure, the further behavior of the 
configuration will be catastrophic. The further evolution of the configuration is not 
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simply that of a live cat, but will be disturbed by two factors: first, by the fact that the 
Gaussian collapse factor will change the shape of | alive), and second, by the fact that 
tails of | dead), which reach the support of | alive) under the Schrodinger evolution, will 
dominate over the contribution from | alive). 

What this example illustrates: First and foremost, this example illustrates how the 
same wave function — the GRW wave function — can be combined with a different PO 
than usual, and thus helps us to get used to the distinction between the wave function 
and the PO. Second, it illustrates how the particle ontology can be combined with dif- 
ferent laws for the wave function — the GRW law instead of the Schrodinger equation 
as in Bohmian mechanics. Third, the example illustrates what it means to derive pre- 
dictions from the PO rather than from the wave function, as the wave function is well 
behaved but the particle configuration is not. In particular, GRWpi forces us to face 
the question: Do the predictions follow from the wave function or from the PO? Fourth 
and finally, GRWpi shows that the GRW wave function can be part of a theory making 
completely different predictions than GRWm and GRWf. 

Another observation: GRWpi is presumably an example of a theory without a for- 
malism. Note first that the \ip\ 2 distribution is not equivariant in GRWpi, in fact there is 
no equivariant density formula at all. This undercuts the reasons for assuming the initial 
configuration was \tp\ 2 distributed, so that we lose the basis for deriving predictions at 
all. But even if we postulated the \ip\ 2 distribution at some point in time (e.g., the big 
bang), so that the theory would make unambiguous predictions, the distribution of the 
outcome Z of an experiment will not be given by a POVM on J4^ ys , presumably not even 
approximately, and we see no reason why the empirical contents could be summarized 
by a formalism at all. 

4.2.2 Bohm's Law and a Modified GRW Law 

To improve GRWpi, one may think of modifying it a bit: Instead of choosing the collapse 
center X at random, as prescribed by the GRW process, one could take 

X = Qj(T) , (32) 

so that the collapse is centered at the actual position of the corresponding particle at 
the time of the collapse. In other words, for every collapse the time T and the label / 
are chosen at random as in the GRW process, but the position of the (center of the) 
collapse is not, but is taken from the particle configuration instead. 

Let us call this model GRWp2; it was called GRWp in [3]. The behavior of a GRWp2 
world is less catastrophic than that of a GRWpi world, but still the \ip\ 2 distribution 
is not equivariant, and so probabilities cannot be expected to agree with \ip\ 2 . As with 
GRWpi, one could consider GRWp 2 with \ip\ 2 distribution at the big bang, but as with 
GRWpi, one presumably obtains no POVM, and no formalism. 

What this example illustrates: The particle ontology can be combined with a mul- 
titude of possible laws for the wave function, each of which is simple and respects the 
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symmetries of the GRW process (invariance under rotations, translations, time transla- 
tions, and Galilean boosts). 

With a little modification in its defining equations, GRWp 2 becomes a better behaved 
theory GRWp 3 [TUl HE]: Instead of (}3"2"j) . take the collapse center X to be 

X = Qj(T) + Z, (33) 

where Z is a random 3-vector that is chosen independently of the past with a Gaussian 
distribution with mean and covariance matrix diag(<r 2 , a 2 , a 2 ). It then follows [48] 
that the conditional distribution of Q(t), given the X, I, T for all collapses up to time t 
(or given ip t ), equals \ipt\ 2 ', that the joint distribution of the ip t for all t > is the same 
as for the GRW process; and that this theory is empirically equivalent to GRWm and 
GRWf. 

What this example illustrates: The empirical content (i.e., the sum of the empirical 
predictions) of the GRWf and GRWm theories can as well be obtained with the particle 
ontology, not only with the flash and matter density ontologies. 

4.2.3 Trajectories Prom the GRW Wave Function 

For the next theory, GRWp4, let us return to the GRW law for ip t , and consider a particle 
ontology with positions given by the wave function by means of the law 

Qi{t) = (i>t\Qi\i>t) , (34) 

where Qi is the position operator of particle number i. That is, the actual position Qi(t) 
is what would in orthodox quantum mechanics be the average position of particle i (if 
measured). In contrast to GRWpi and GRWp 2 (and Bohmian mechanics), this theory 
does not require any initial data about the particle configuration, as the configuration 
is a function of ip t - Like GRWpi and GRWp 2 , this theory GRWp 4 is not empirically 
adequate. 

This is best seen when considering a system of identical particles, whose wave func- 
tion obeys permutation symmetry or anti-symmetry. As a preparation, we note that the 
GRW process as we have defined it in Section 12.11 above does not preserve permutation 
symmetry. For identical particles the GRW process needs to be defined differently to 
include a collapse mechanism which preserves permutation symmetry; for the technical 
details (which will not matter for our purposes) see [44]. So consider the theory defined 
by that process together with (134"]) . For a wave function with appropriate permutation 
symmetry, fl34"]) implies that all particles (of the same species) have the same position. 
However, it is empirically incorrect that all electrons have the same position. 

What this example illustrates: Apart from being an example of the multitude of 
possible laws for the PO, an interesting trait of this theory is that the configuration of 
the PO supervenes on the wave function by means of the law (1341) . as it does in GRWm 
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and GRWf. GRWp4 (or the version of it for identical particles) shows that among the 
many ways in which a configuration of matter (be it a particle configuration, a continuous 
matter distribution, or flashes) can supervene on the wave function, different possibilities 
may strongly disagree about the empirical predictions. 



4.2.4 Configuration Jumps and GRW Law 

Another theory, GRWps, has the following laws: The wave function ip t follows the GRW 
process, and the configuration moves according to Bohm's law of motion (T29|) between 
the GRW collapses. However, at the time T when ip collapses around X e R 3 with label 
J, also the configuration Q jumps; more precisely, only the 7-th particle jumps, and it 
jumps to the random center X of the GRW collapse: 

Qj(T+) = X. (35) 

Again, \ip\ 2 is not equivariant, and indeed, the behavior of the particles is quite 
catastrophic: If ip(T—) is the wave function of Schrodinger's cat, the configuration 
Q(T—) is that of a live cat, and the wave function collapses to that of a dead cat, it 
may do so with just a few collapses connected to a few particle labels, corresponding to 
particles that would have to be in different positions depending on whether the cat is 
dead or alive. As a consequence, the configuration after these few collapses will be one 
of a live cat, with a few particles moved to where they would have to be if the cat were 
dead. So, this configuration is very different from what one would normally associate 
with | dead). Also, this configuration may be well outside the support of both |dead) and 
| alive); it may be a configuration for which \ip\ 2 is literally zero, or much smaller than 
even the remains of | alive) after the collapse. And the behavior of such configurations 
should be expected to be catastrophic. 

What this example illustrates: One normally thinks that when ip = |dead) then there 
is a dead cat. However, as GRWp 5 illustrates, in the PO view this cannot be taken for 
granted but must be checked. 

4.2.5 Another Way of Configuration Jumps and GRW Law 

Within the framework that tpt follows the GRW process and the configuration moves 
according to Bohm's law of motion f l29|) between the GRW collapses, further options 
besides GRWps come to mind. Instead of f )35|) . we may postulate that at the time T 
when ip collapses around X e M 3 with label /, the 7-th particle jumps to a random 
position distributed with density 

P(Q/(T+) € dx) = \^t+ (gig) • • • Qj-i(T), x, Q I+1 (T) . . . Q N (T))\ 2 

dx ~ J dx'\^ T+ (Qi{T) . . . Q/_i(T), x', Q I+1 (T) . . . Q N (T))\ 2 ' ° '* 

which is the I'f/'T+I 2 distribution conditionalized on the configuration of the other parti- 
cles. 



Another possibility, which we call GRWp6, is that not just the 7-th particle (the one 
associated with the collapse) jumps, but all particles jump. Specifically, choose Q(T+) 
at random with distribution |^ + | 2 . 

In GRWp6, the |?/>| 2 distribution is indeed equivariant, in the sense that the con- 
figuration Q(t) will always have distribution \ipt\ 2 - As a consequence, this theory is 
presumably empirically equivalent to GRWm and GRWf (in the sense that there is no 
experiment that could distinguish GRWp6 from GRWm and GRWf). However, the par- 
ticles in GRWp6 do not necessarily behave in a reasonable way: For example, consider 
an agglomerate of N = 10 11 particles and arrange a superposition of two well-separated 
locations, ip = ■^(|here)® Ar + Ithere)®^). Suppose the particles are all "here" (this 
happens with probability 1/2) and that the first collapse, which occurs after about 1 
day, reduces if) to |there) 0Ar (this happens independently with probability 1/2). In this 
case all particles jump from "here" to "there." 

What this example illustrates: Even if a theory is empirically equivalent to a reason- 
able theory (such as GRWm and GRWf) it need not itself be a reasonable theory]^ 

4.3 MBM: Bohm-Like Trajectories From the Master Equation 

The most interesting example in our list of toy theories is perhaps MBM; the abbre- 
viation stands for "master equation Bohmian mechanics." MBM resembles Bohmian 
mechanics in that it is deterministic and that its PO consists of particles, but at the 
same time it is empirically equivalent to GRWf and GRWm, as we will show in Sec- 
tion gxu 

The law of motion (|29|) is replaced by the following equation (considered already in 
[26] ) using, in the role of the wave function, a density matrix p: 

d Qk P(n n x h T V qk (q\p\q') 

~ %{Qi,...,Qn) = Im . (37) 



dt *v> mk {q \ p \ q >) 



9=9'=(Ql ,->Qjv) 



The density matrix p evolves according to (EJ), which we repeat here for convenience: 

^ = -i[H,p t ] + \J2 I d 3 xgt /2 x p t gl /2 x -NXp t . © 
k=i J 

We make a few comments on how these equations are to be understood. 

Eq. (I37p is the natural generalization of Bohmian mechanics to density matrices, and 
reduces to (I2"9~j) in case of a pure state p = However, it is important to notice 

that the density matrix considered here is not the one that is normally regarded as the 
density matrix of a system, which arises by averaging \if))(ip\ (i n case the wave function 
if) is random) or by tracing out the environment of the system. The density matrix in 



4 This point is also illustrated by the following theory [14] that is empirically equivalent to Bohmian 
mechanics: Using a Schrodinger (i.e., non-collapsing) wave function tp t and a particle ontology, let the 
configuration Q(t) be random with distribution \ip\ 2 , independently of the past. 
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fl37|) . in contrast, does not arise from averaging or partial traces but is assumed to be 
one of the fundamental variables of the theory. The complete description of the state 
is, instead of the pair (Q,i/j) in Bohmian mechanics, the pair (Q,p). 

The master equation ([7]) is, in the GRW theories, a consequence of the GRW evo- 
lution of the wave function. This is different in MBM. In MBM there is no random 
wave function. In MBM, ([7]) holds by fiat, not as a theorem. The defining equations of 
MBM— its postulates— are (17j) and (13Tj). 

What this example illustrates: MBM shows that the empirical content of GRWm and 
GRWf is compatible with a deterministic theory, and in particular does not imply wave 
function collapse: after all, MBM involves the master equation ((TJ but not literal wave 
function collapse as in ([1]). If experiments will confirm the GRW deviations from quan- 
tum mechanics, then Bohmian mechanics can be modified so as to reproduce the GRW 
predictions. MBM also illustrates how it can make sense to speak of the density matrix 
of the entire universe and, more specifically, how a density matrix can be a fundamental 
object and part of the ontology, rather than just encoding statistical information. 

4.3.1 Empirical Equivalence of MBM with GRWm and GRWf 

We now prove the empirical equivalence. The derivation of predictions from MBM is 
analogous to that from Bohmian mechanics. The analogue in MBM of the quantum 
equilibrium distribution p^ described in ( 13T|) is the distribution 

P P (q) = (q\p\q) • (38) 

Note that p p (q) > 0, and fQP p {q)dq = tip = 1. As we show in Appendix IA~1 this 
distribution is equivariant in MBM as a consequence of fl37|) and ([7]). It is therefore 
consistent to assume, as we do, that the configuration Q t of the world has distribution 
p Pt at every time t. And therefore, the probability at time t of a certain macroscopic 
configuration is 

p(S)= [ dq(q\p t \q)=ti(p t P(S)) (39) 
Js 

where S is the set of all microscopic configurations consistent with that macroscopic 
configuration, and P(S) the projection operator corresponding to the set S, defined by 

P(S)= [ dq\q)(q\. (40) 
Js 

In GRW theories, p(S) can be written in terms of the GRW wave function ip t as 

tr(pt P(S)) = [ n^t E dcf>) \\P(S)<P\\ 2 . (41) 

Since ip t is typically concentrated on a single macro- configuration, the probability dis- 
tribution P(ipt £ dip) is typically concentrated on those <fi with either ||P(S')0|| ~ or 
||P(S')0|| ~ 1; thus, p(S) equals the probability that tj) t is (nearly) concentrated on S. 
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And in this case, either the flashes of GRWf or the matter density of GRWm gives rise 
to the same macroscopic appearance as configurations from S. In other words, at any 
fixed time the MBM, GRWf, and GRWm worlds have the same probability distribution 
over the possible macro-states. 

Now empirical equivalence follows immediately: If there were an experiment which 
had (probably) one outcome Z = z\ in MBM and another one, Z = z 2 m GRWm and 
GRWf, then at the time when the experiment is finished, the probability distribution 
over the macro-states would have to be different in MBM than in GRWm and GRWf, 
but it is not. 

As a consequence of the empirical equivalence between MBM, GRWm, and GRWf, 
the empirical content of MBM is summarized by the GRW formalism [31]. (In fact, the 
GRW formalism was first discovered starting from MBM.) Our argument concerning 
empirical equivalence also exemplifies that empirical equivalence is a statement about 
the PO, as discussed in Section [3j 

4.4 Master Equation and Matter Density 

If one can consider a version of Bohmian mechanics in which the density matrix plays 
exactly the role of the wave function, then why not do the same trick with the matter 
density and flash ontologies? 

For the matter density ontology, this would mean to postulate, in analogy to and 
replacing ffTOl) . 

m(x, t) = tr (p t M(x)) (42) 

with M(x) = ^2 i m i 5(Q i — x) the mass density operator, as before. Here, p t is taken 
to evolve according to the master equation We thus obtain a theory that could 
be called Mm, in which the fundamental objects are a density matrix (which, as in 
MBM, does not represent an ensemble, or the observer's limited knowledge, but is, by 
postulate, a fundamental object) and the continuous matter with density m(x,t). This 
theory, though, is very different from GRWm! It has a many-worlds character. For 
example, if at some initial time p = \tp)(tp\ with = (|dead) + | alive) )/\/2 being the 
wave function of Schrodinger's cat, then after a short while the GRW function ip t will be 
either |dead) or | alive), but p t will be | dead) (dead | + | alive) (alive |, up to a factor |. As a 
consequence, the m field of GRWm will be either mdead or m a i ive , but the m field of Mm 
will be m dead + m alive , up to a factor ~. Both cats are there at once, but with reduced 
mass (which the cats, however, do not notice). A very similar theory "Sm", with the 
unitary Schrodinger evolution instead of the master equation ([7]), has been described 
in some detail in [31 0]. It seems possible, perhaps even likely, that Sm is empirically 
equivalent to standard quantum mechanics [3|; if so, then for the same reasons Mm 
should be empirically equivalent to GRWm and GRWf. 

What this example illustrates: Foremost, this example illustrates the big difference, 
for quantum theories without observers, between a density matrix and a random wave 
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function. Since with every probability distribution over wave functions there is asso- 
ciated a density matrix, and since for the purpose of computing predictions only the 
density matrix is relevant, it is common practice in quantum mechanics to immediately 
replace every probability distribution over wave functions by the density matrix. Here, 
however, the distinction is crucial: The deterministic m function obtained from the 
density matrix provides a many-worlds picture of reality, whereas the random m func- 
tion obtained from a random wave function that is either |dead) or | alive) provides a 
single- world picture. 

4.5 Master Equation and Flashes 

What happens to GRWf when we replace the wave function by a density matrix? Two 
versions of what this might mean come to mind: 

MGRWf: Postulate that the every time i, a density matrix p t , whose 

evolution will be described below. Postulate further that flashes occur, as in GRWf, at 
random times T with constant rate iVA, and for a random I £ {1, ... , N}. In contrast 
to GRWf, the probability density of the flash location X is given by 

P(x) = -^—^ = tr (p T - gi, x ) (43) 

instead of \\ipT- 9i1 1| 2 in GRWf. Postulate further that when a flash occurs at time T, 
the density matrix changes to 

1 V 2 V 2 (AA\ 

Pt+ = -^glx PT- 9l' X ( 44 ) 

with C a normalizing constant. Between collapses, pt evolves unitarily as usual. Note 
that if p = | "00 ) ("00 1 is a pure state then it will remain pure and in fact evolve according 
to the usual GRW dynamics. 

It follows that the joint distribution of all flashes is given by the probability measure 

PmgW^ eS)= tr (p G (S)) (45) 

with S any set of flash histories and Gq(-) the POVM governing the distribution of the 
flashes in GRWf, see fl2Tj) . As a consequence of (Jl5j) . the theory MGRWf, although 
formulated in terms of a density matrix, is physically equivalent to GRWf! 

Before justifying this claim, let us elucidate the notion of "physical equivalence," 
which can be given a clear definition in the framework of PO [3]: Consider a theory 
T (e.g., GRWf with a particular choice of A, a, N, and potential V) and initial data 
.Do for T (e.g., for GRWf, an initial wave function ip ); together, T and D define a 
probability distribution P over histories of the PO after the initial time (e.g., for GRWf, 
a distribution over flash patterns in space-time). Now we say that the pair (T',D' Q ) is 
physically equivalent to (T, D ) if it defines the same distribution P as (T, D ). That is, 
two descriptions of a universe are physically equivalent if they provide the same history 
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of the PO, or (if appropriate) the same probability distribution thereof; variables (such 
as ip) that are not part of the PO may be different in the two descriptions. 

For the purpose of comparison with MGRWf, suppose that, as the initial data D for 
GRWf, we do not specify the initial wave function but instead a probability distribution 
fi over initial wave functions. That still defines a distribution P over the PO histories: 
Since the joint distribution of all flashes, given ip , is (■0o|G ? o(-)|V ; o), the unconditional 
joint distribution of all flashes is 

^GRWf(FeS)= J »(di/>)(1> \G (S)\il>o)=tr(p G (S)), (46) 

with p the density matrix of the ensemble p. Since this is the same formula as (|45p . 
we obtain the physical equivalence. That is, the further possibilities of initial data p in 
MGRWf do not lead to new flash histories. 

How is MGRWf related to the master equation ([7])? Since p t is random, it does not 
evolve according to (J7J). However, there is another density matrix, namely 

Pt = Ept , (47) 

where E means the expectation over the random flashes; p t does evolve deterministically 
according to the master equation (j7J. 

What this example illustrates: First, it illustrates the concept of physical equivalence: 
There is no physical difference between GRWf with a random initial wave function of the 
universe and MGRWf, as the distribution of the PO is the same. Second, it illustrates 
the different roles that a density matrix can play: On the one hand, it can be part 
of the ontology as one of the fundamental objects, such as the random density matrix 
p t . On the other hand, the other density matrix p t is a mathematical object encoding 
information about the probability distribution of pt- Thus, finally, it also illustrates 
that it is not necessarily appropriate to speak of "the" density matrix. 

Mf: Another way of replacing the wave function in GRWf by a density matrix is to 
retain (|43p for the distribution of the flash location X but not adopt the collapse rule 
( 144"|) . Instead, let p t evolve deterministically according to the master equation (J7J. In 
other words, Mf arises from MGRWf by replacing p t in ( 143]) by pt. 

This theory, too, has a many-worlds character, as the set F of flashes will be the 
union F = F deild U -F alive of a set of flashes of a live cat and a set of flashes of a dead cat 
(similar to the model "Sf" considered in [4j). For the same reasons as for Mm, Mf is 
presumably empirically equivalent to GRWf and GRWm. 

What this example illustrates: By way of contrast with MGRWf, it illustrates the 
difference between the "collapsing" density matrix of MGRWf and the deterministic one 
arising from the master equation (J7J). It thus illustrates why it is important that the 
law governing the PO be precisely formulated (since in particular, it specifies precisely 
which density matrix to use). 
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5 Conclusions 



In this paper we have illustrated the notion of PO (primitive ontology, i.e., variables 
describing the distribution of matter in space and time) and its usefulness for a clean 
derivation of empirical predictions. We have done so by (i) describing several toy exam- 
ples of theories with a PO (Section HJ), (ii) studying their predictions (Section Hj), (iii) 
reviewing known theories with a PO (Sections 121 and l4.ip . and (iv) reviewing known 
derivations of predictions (Section [3]). In so doing we have illustrated the different roles 
that the PO, the wave function, and the density matrix have in such theories. 



A Proof of Equivariance in MBM 

We show that equivariance of p p (q) = (q\p\q) follows from (j7|) for the Hamiltonian (j3J) 
and fl37j) . Before we give the formal proof, we note that the essential reason is that the 
non-unitary (diffusive) terms in (J7J) (i.e., the second and third term on the right hand 
side), do not contribute to the continuity equation for p p . This can be understood by 
noting that a GRW collapse does not change the diagonal entries (q\pt\q) of the density 
matrix in the position representation. 
Here is the equivariance proof. 
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Note that the potential V in H does not contribute in (149!) because 



<?I^|?> 



/ (q\V\q') (q'\pt\q) = V(q)(q\p t \q) = (q\ptV\q) . 

=S(q-q')V(q) 



(55) 
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Since any probability distribution p on configuration space will be transported, under 
the flow (1571) . according to the continuity equation 

% -EV ft 'K)> (56) 

k=l 

we have that p t = p pt satisfies ( 1551) . which is equivariance. 
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